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ABSTRACT 

We present results from three-dimensional (3-D) numerical simulations investigating 
the magnetohydrodynamics of cloud-wind interactions. The initial cloud is spherical 
while the magnetic field is uniform and transverse to the cloud motion. A simplified 
analytical model that describes the magnetic energy evolution in front of the cloud is 
developed and compared with simulation results. In addition, it is found the interaction 
of the cloud with a magnetized interstellar medium (ISM) results in the formation of 
a highly structured magnetotail. The magnetic flux in the wake of the cloud organizes 
into flux ropes and a reconnection, current sheet is developed, as field lines of opposite 
polarity are brought close together near the symmetry axis. At the same time, magnetic 
pressure is strongly enhanced at the leading edge of the cloud from the stretching of the 
field lines that occurs there. This has an important dynamical effect on the subsequent 
evolution of the cloud, since some unstable modes tend to be strongly enhanced. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - magnetic fields 



1. Introduction 

Magnetic fields are a pervasive element of the interstellar and intergalactic medium and they 
are often extremely relevant in characterizing local and global behaviors of these media. Mag- 
netohydrodynamic (MHD) cloud-wind interactions are believed to be important in determing the 
observed filamentary and clumpy morphology associated with clouds moving through a magnetized 
interstellar medium. These processes may also result in a local enhancement of the background 
magnetic field, which, in turn, provides an important feedback on the subsequent evolution of the 
cloud. Jones et al. (1994, 1996) have shown that these regions of strong magnetic pressure (what we 
call "magnetic bumpers" ) develop in front of the cloud as a result of the stretching of the field lines 
that anchors on the cloud surface. Such regions may also be adjacent to shocks, which can serve 
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as acceleration sites of high energy particles (e.g., Jones & Kang 1993). Within some supernova 
remmants there is also clear evidence for supersonic clumps (e.g., Jones et al. 1998). Consequently 
the strong field regions can show enhanced nonthermal radio emission (Jones et al. 1994). Miniati 
et al. (1997, 1999a) investigated the role played by magnetic fields in cloud collisions by comparing 
two dimensional (2-D) hydrodynamic and magnetohydro dynamic simulations. They concluded that 
a magnetic field transverse to the cloud motion can dramatically alter the outcome of the collisions, 
preventing the disruption of the clouds otherwise occurring in almost all other scenarios. Since the 
development of the magnetic bumper was crucial for this result, Miniati et al. (1999b) have extended 
the Jones et al. (1994, 1996) work. In particular, they further investigated the formation of such a 
bumper under a broader range of initial conditions by studying the propagation of clouds through 
various oblique magnetic fields. In addition they also assessed the issue of the exchange of magnetic 
and kinetic energy during the evolution of diffuse clouds in the ISM. There, in fact, they offer a 
possible explanation for the comparable values of magnetic and kinetic energy densities observed in 
some H I complexes (Heiles 1989; Verschuur 1989; Myers & Khersonsky 1995; Myers et al. 1995). 
Gregori et al. (1999) have presented three dimensional (3-D) numerical simulations that show 
dramatic dynamical effects of the magnetic field in determining the cloud evolution during its 
propagation through a magnetized medium. In particular, those authors found that a strong field 
enhances the development of Rayleigh- Taylor unstable modes, thus hastening the cloud disruption. 
On one hand this demonstrates that magnetic fields in 3-D do not simply slip around the cloud 
surface, somewhat similar to the 2-D case. On the other hand, the 3-D influence of the field was 
disruptive, opposite to that seen in 2-D simulations. 

In this paper we present 3-D numerical simulations of moderately supersonic cloud motion 
in a magnetized interstellar medium. The cloud is treated as non self-gravitating and adiabatic. 
First we present a detailed analytical model for the stretching mechanism of the field lines and the 
consequent magnetic field amplification at the cloud nose. We also show that the final outcome in 
a cloud-wind interaction is the development of complex features analogous to the ones observed in 
cometary plasma tails, resulting in the formation of a highly structured magnetotail. The magnetic 
field itself organizes in coherent tails, or flux ropes, in the wake of the cloud. These ropes are 
associated with the development of a Sweet-Parker reconnection sheet that alters the field topology 
there. Our aim is to give a picture of the basic processes that develop in cloud-wind interactions, 
in order to provide useful insight for observations in a large variety of astrophysical environments 
(see e.g., Dgani & Soker 1998). The use of 3-D simulations is clearly a considerable advantage 
over previous work, since now the full spatial domain can be investigated without the geometrical 
limitations imposed by 1-D or 2-D calculations. 

The paper is organized as follows. In §2 we describe the numerical setup, and the character- 
istic physical parameters of the problem. The computational results are introduced in §3. The 
development of the Rayleigh- Taylor instability and the cloud disruption are briefly discussed in 
§3.1. The analysis of the magnetic energy evolution is given in §3.2. In section §3.3 we investigate 
the formation of flux ropes in the wake of the cloud and in §3.4 the magnetic reconnection. Our 
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results are summarized in §4. 



2. Numerical Setup and Definition of the Problem 

The numerical computation is based on a total variation diminishing (TVD) scheme for ideal 
MHD (Ryu & Jones 1995). This is an explicit, conservative finite-difference method with second 
order accuracy in space and time. We have used the multidimensional, Cartesian version of the 
code (Ryu, Jones & Frank 1995) with a constrained transport scheme for preserving V • B = 
(Ryu et al. 1998). Neglecting self-gravity and radiative energy losses, the complete set of the 
simulated equations for the velocity (u), magnetic field (B), density (p), pressure (p) and "color 
tracer" (C) can be conveniently written as 

^ + V • (pu) = 0, (2-1) 
dw 1 1 

v + u-Vu + -Vp--(VxB)xB = 0, (2-2) 
at p p 

^ + u • Vp + 7pV • u = 0, (2-3) 
i9B 

— - V x (u x B) = 0. (2-4) 
at 

dC 

— + u • VC = 0. (2-5) 
at 

The last equation has been added to the standard set of ideal MHD equations in order to be able 
to follow the motion of the cloud material itself. C corresponds, in fact, to the mass fraction of the 
cloud gas inside the computational cell. 

Initially, all the cloud material is labeled with C = 1, and the ambient medium with C = 0. 
At any time t, the density of cloud material in a fluid cell is given by p c = pC, where p is the 
total fluid density at that point. In the MHD equations, the magnetic field is normalized such that 
the factor Att does not appear, giving an Alfven speed va = B/y/p. We assume an adiabatic index 
7 = 5/3, initial pressure equilibrium at po = 3/5, and an initial density in the background medium 
p = pi = 1- Thus, the velocity is expressed in units of the sound speed in the ambient medium: 
c s = ilPo/ Pi) 1 ^ 2 = 1- The initial cloud density is p c = xPi> with \ = 100- A thin transition layer 
~ 0.2R C (R c is the cloud radius) around the cloud, introduced to reduce the Richtmyer-Meshkov 
instability at startup, brings the density to the value of the intercloud medium. A finite transition 
layer is in general expected due to both thermal condution (Balbus 1986) and photoionzation 
(Tielens & Hollenbach 1985) at the cloud boundary, although at this point in our simulations it is 
considered only for reasons of numerical stabililty. The cloud is initially spherical in shape. The 
numerical value for its radius, R c , is set to unity, and this is chosen as the unit of length. This also 
sets the unit time to the cloud sound crossing time, r cs = R c /c s (= 1 in numerical units). At time 
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t = 0, the cloud is set in motion with respect to the uniform background medium. Its velocity is 
u c = Mc s , with a starting value for the intercloud Mach number, M = 1.5. The magnetic field is 
conveniently expressed in terms of the familiar parameter 



where ps = B 2 /2 is the magnetic pressure. In our numerical simulations we have considered both 
the cases of an initially strong field ((3 = 4) and a weak field ((3 = 100). To be able to compare 
with pure hydrodynamic effects, a case with (3 = oo (no magnetic field) has also been computed. 
A summary of all the simulations performed is given in Table 1. In addition to the (3 parameter, 
another important dimensionless number often used to describe the action of the magnetic field on 
the fluid motion is the Alfvenic Mach number Ma = u/vai, where vai = B / \fpi is Alfven velocity 




in the intercloud medium. So, we initially have Ma = 2.74 for to (3 = 4, and Ma = 13.7 for 



The computational domain is outlined in Fig. 1. Symmetrical boundary conditions are em- 
ployed on the y = and z = planes, inflow conditions are applied on the x = plane, while open 
conditions are used on all other boundaries. This choice eliminates odd modes of instabilities, but in 
companion, fully 3-D simulations at somewhat lower resolution we saw no evidence that such modes 
play a deciding role in cloud evolution. We employed a uniform grid, N x x N y x N z = 416 x 208 x 416, 
spanning i the volume of interest. The volume computed was bounded along x-y, x-z planes 
through the initial cloud center. This gives a resolution of 26 zones per cloud radius, less than 
that in our previous 2-D MHD simulations (Miniati et al. 1999b, Jones et al. 1996). So, small scale 
surface perturbations were relatively more damped. However, from our experience, the adopted res- 
olution is sufficient to capture basic cloud evolution over the time interval considered. In addition, 
we have carried out several lower resolution 3D simulations spanning the full volume of interest 
(see Table 1). Some aspects of those companion simulations are reported as well in the following 
sections. The cloud shapes (sphere and cylinder) that we have considered are quite ideal and they 
need to be considered to give only a qualitative picture of the true interaction of a magnetized 
wind with an interstellar cloud which has no distinctive shape. In this respect, the simulations 
with spherical and cylindrical clouds show a quite similar behavior. Tests carried out with a dif- 
ferent shape, an elliptical cloud, have also confirmed the same pattern. They all show qualitative 
agreement with the simulations reported here. 

Since the cloud motion is supersonic, its motion leads to the formation of a forward, bow shock 
and a reverse, crushing shock propagating through the cloud. The approximate time for the latter 
to cross the cloud is referred to as the "crushing time" 4 (e.g., Jones et al. 1994): 




(2-6) 




= 100. 



2R cX 1/2 



(2-7) 



Mc s 



4 This form of the crushing time is a factor of 2 larger than the one of Klein et al. (1994), since our definition is 
based on the cloud diameter instead of the cloud radius. We use this definition since it more closely measures the 
actual time before the crushing shock emerges. 
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Since the crushing time corresponds to the typical scale for the cloud evolution, in the following 
figures and discussions, time will be expressed in terms of r cr . 

Finally, it is necessary to add a few words about the significance of our choice in the initial 
direction of the field lines. In our simulations, the initial magnetic field has been set up transverse 
(perpendicular) to the cloud motion. As pointed out in previous 2-D simulations (Miniati et 
al. 1999b; Jones et al. 1996; and Mac Low et al. 1994), a magnetic field aligned with the direction of 
the cloud motion never becomes directly dynamically relevant in terms of body forces, even if it may 
have some stabilizing effects. For a general orientation, there will always be some component of the 
field transverse to the cloud velocity, which will be stretched around the cloud body. In this respect, 
our simulations can be viewed as an approximate solution of the more general problem of the oblique 
field orientation. As stressed by Jones et al. (1996) and Miniati et al. (1999b), for supersonic bullets 
most field directions will indeed produce effects similar to the transverse field case. To confirm that 
the field evolution in the experiments described here are not special cases geometrically, we have 
also carried out fully 3-D simulations at low resolution, including oblique initial field orientations, 
confirming the general behaviors for the more restricted symmetries imposed. 



3. Results 

3.1. Dynamical Evolution & Cloud Disruption 

As discussed in a companion paper (Gregori et al. 1999), the magnetic field amplication that 
occurs in front of the cloud inhibits instabilities in the x-y plane, but hastens those in the x-z plane, 
thus accelerating the process of cloud disruption. An example of the effects of this magnetically 
enhanced Rayleigh- Taylor (R-T) instability can be seen in Fig. 2. This shows density slices in the 
y = plane for the (3 = oo (top), (3 = 100 (middle) and (3 = 4 (bottom) simulations, at t = 0.94r cr 
(left) and t = 2.25r cr (right). Clearly, for the (3 = 4 case, the growth of a R-T instability is much 
more rapid and, unlike the other cases, very pronounced density fingers have developed at the 
simulation end. In addition, to illustrate how the cloud shape evolves, we define generalized cloud 
sizes in terms of the coordinate moments of inertia. Following Klein et al. (1994) and Xu & Stone 
(1995), the cloud extension in the ith direction, at time t is given by 

where r% is the ith position coordinate with respect to the center of mass, and the integrals are 
intended over the entire computational domain. In Fig. 3 we have plotted R y and R z for both 
the MHD and the hydrodynamic simulations. Compared to the hydrodynamic case, the wrapping 
of the field lines around the cloud in the x-y plane produces a strong radial magnetic pressure 
gradientgradient that squeezes the cloud gas. Such pressure, consequently, forces an extrusion of 
the cloud along the z direction. This explains the trends visible in Fig. 3. There we can read that 
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for the (3 = 4 case, after one crushing time, R y has already been reduced by 30-40% while R z has 
been increased by a similar amount. We may observe that at t > 1.5r cr , the expansion in the z 
direction is sustained at a very large rate, while R y stays almost constant. Finally, in the (3 = 100 
simulation the change in the cloud form follows the same qualitative pattern as for the (3 = 4 case, 
although the evolution is considerably less dramatic and rapid. 

3.2. Magnetic Energy Evolution 

3.2.1. Model 

In the attempt to understand the process of the magnetic bumper formation and its inter- 
action with the surrounding flow, Miniati et al. (1999b) and Jones et al. (1996) have proposed 
an approximate model based on the Faraday induction equation. Their analysis showed that the 
magnetic energy first increases with time exponentially (Jones et al. 1996) and then according to 
a power-law of index between one and two (Miniati et al. 1999b). In hopes that it may serve as a 
simple tool for understanding field growth in such problems, we examine again the magnetic energy 
evolution explicitly using the flux freezing condition while incorporating the important and more 
realistic assumption of a limited region of velocity shear around the cloud nose. 

In fact, the time evolution of the magnetic field, in ideal MHD is given by the solution of the 
following equation in Lagrangian form 

dt ( p ) p VU 
This equation can be formally integrated as {e.g., Batchelor 1967) to give 

(!)-(f).-^ 

Here, (B/p)o refers to the value at time t = 0, and a is the initial position vector of the fluid 
element being followed. It is worth stressing that the previous equations are only valid far from 
shocks. This is indeed the case in our analysis since the magnetic energy tends to evolve mostly in 
front of the cloud, far from the bow shock and the crushing shock penetrating into the cloud. The 
mapping function of the flow field is X(a, t), which locates the fluid element at subsequent times. 
Its derivative, 9X/<9a, parametrically measures the deformation of a fluid element, and, therefore 
determines the local growth of the magnetic field. From eq. (3-3) the strength of the magnetic field, 
|B| = B, then evolves according to: 

fB\ (BA d v . . (B\ dX . 
where X = |X| and a = |a|. In order to derive the evolution of the term dX/da we assume steady 



(3-2) 



(3-3) 
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flow and write the equation for a fluid element along a flow line as (e.g., Aris 1962): 

dX ~ f ^ + ^ cos sin dX, (3-5) 



d(dX) dui dxj dxi f du x du y 



dt dxj dX dX \ dy dx 

where we neglect the terms involving u z and indentify dxj/dX and dxi/dX as the direction cosines 
of the flow line. To integrate equation (3-5) we make the important and realistic assumption 
that the velocity shear vanishes outside a finite region in front of the cloud. This assumption is 
confirmed by Fig. 10 where we can see that in most cases the velocity field is characterized by 
a strong shear only around the cloud front. We emphasize that if a fixed portion of a fluid line 
is subject to stretching, then its length, and, hence, its strength, will grow only linearly (as shown 
below by eq. 3-6), as opposed to an exponential increase occuring when the stretching is constant 
over the full extent of the fluid line. In addition, we assume for simplicity that such a shear pattern 
is approximately constant in space. Therefore, integrating eq. (3-5) in both space and time we 
obtain in that restricted region: 

v , \ /"* , f a+la (du x du v \ n . n , v 1„ (du x du v \ . . 

where i a represents the length of the flow line over which the velocity shear is non-null and the term 
sin cos 6 has been taken to contribute a factor \ coming from integrating along an arc size of the 
order of ^, as appropriate for the cloud nose. We are now able to derive the quantity of interest, 
namely dX/da. Assuming a laminar flow regime before the cloud and around its nose, it is possible 
to see that the change of the length (along which the strain occurs) S£ a between two adjacent flow 
lines is of the order of their separation 8a, or, equivalently 5£ a /8a = 1. We approximate also the 
velocity shear \ (dv x /dy + dv y /dx) ~ X(t)u c /R c , since the flow speed increases from zero on the 
cloud nose to roughly the cloud speed, u c , along an arc of length ~ R c at the start of the simulation. 
Here, X(t) is a slowly varying quantity introduced to correct for all the details of the MHD flow 
dynamics that have not been included in this simplyfied analysis. From equations (3-4) and (3-6) 
we have then 

B B ( u c \ , N 

- ~ 1 + X-^t , (3-7) 



P PO V R c / 

or in terms of the magnetic energy density (magnetic pressure) ps = B 2 /2, 

™°(£) 2 ( 1+A i ( ) 2 ' (3 - 8) 

The quantity X(t) parameterizes the evolution in the flow. In the beginning this mostly represents an 
expansion of the flow field in response to increased magnetic pressure; that is, the flow lines become 
more spread out. Eventually, however, as the cloud begins to decelerate, it includes reductions in 
its asymptotic speed, as well. The details of the cloud dynamics embodied in X(t) are difficult to 
estimate. Nevertheless, based on the results of our 3-D simulations and on the theoretical model of 
Miniati et al. (1999b), we believe that for t > t s = R c /u c a functional form X(t) ~ (t/r s )~ q (with 
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q ~ 5) is adequate to capture the main features in the magnetic field evolution. Then, we can 
rewrite eq. (3-8) as becoming asymptotically 

i ™(£) 2 (£) 2M (3 - 9) 

For q = I the magnetic pressure would continue to grow linearly with time. As discussed below, 
our simulations indeed suggest that the field pressure tends to increase with a power law exponent 
of order unity, at least until the crushing shock has completely crossed the cloud. 



3.2.2. Quantitative Analysis and Comparison with Numerical Simulations 

Neglecting the initial times (t < r s ), when compressibility is dominant, in the following we 
will assume p ~ po and only consider the time regime in which equation (3-9) is appropriate in the 
flow around the cloud. Then, let us estimate the integral form of equation (3-9) over the entire 
computational volume V. We indicate with V e the region around the cloud where the velocity shear 
is large, and, therefore, field amplification takes place. Then Xu c /R c becomes the average rate of 
strain over such a volume, outside of which the magnetic field remains approximately equal to its 
initial value. Typically, V e is of the order of V c , the cloud volume. 

E B = j^PBdV ~p B0 V c ^-\ + PB0 (V-V C ). (3-10) 
Then, we can rearrange (3-10) as 

Eb — Ebo AEb 



^B0 ^B0 



(3-11) 



where m = 2(1 — q), Ebo = PboV is the initial magnetic energy in the computational volume 
and EgQ = pboVc is the initial cloud magnetic energy. In Fig. 4 (top panel) we plot the time 
dependence of AEb/E^q obtained directly from the numerical simulation to compare with this 
relation. A lower resolution simulation result (16 zones per cloud radius) is also given for the 
(3 = 4 case. Both the high and the low resolution curves are sufficiently close to suggest a good 
convergence for the global magnetic energy evolution. It is clear that the total magnetic energy in 
the volume tends to increase monotonically over the simulated time. 

Fitting the numerical results in the time interval 0.01 < t/r cr < 1.0 to a t m dependence, gives 
a power law index m pa 1.1 for [3 = 4 and m pa 1.4 for [3 = 100. These results suggest m — ► 1, at 
least during this interval in the cloud evolution. Indeed, a linear increase in the magnetic energy is 
to be expected. In fact, up to t ~ r cr the cloud does not suffer any deceleration. Since the system 
is then supplied with a steady flux of kinetic energy, we may expect that a constant fraction of the 
latter is converted into magnetic energy at a constant rate. However, Fig. 4 shows that, in both 
cases, at t > r cr the rate of magnetic energy growth drops below linear. This is a signature of the 
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beginning of the cloud deceleration and the consequent reduction in energy supply, although some 
outflow of magnetic energy from the computational volume is also occuring (see below). It is of 
interest to compare the increment in magnetic energy to the evolution of the total kinetic energy 
inside the computational box. This is also shown in Fig. 4 (bottom panel). From it we learn that 
at about one crushing time more than 10% of the initial kinetic energy has been converted into 
magnetic energy for the (3 = 4 case. This fraction becomes larger that 15% at simulation end. On 
the other hand, for the (3 = 100 case, the conversion amounts to only a few percent, in accordance 
with a much larger conversion timescale as determined later (see §3.2.3). We point out that the 
increments in the magnetic energy correspond to and are responsibile for the decrements over time 
of the kinetic energy (measured with respect to an observer who sees the intercloud medium at 
rest) for both values of (3 as plotted in Fig. 5. In reality, a closer analysis reveals that at simulation 
end, the decrements in kinetic energy are somewhat larger than the correspondent magnetic energy 
increments, indicating that some magnetic energy flux has escaped the computational box. Thus, 
the final values for the AEb in Fig. 4 are understood only as lower limits. Comparing results 
at different resolutions for the (3 = 4 case, again indicates a good convergence of the simulations 
in this property. It is worth mentioning that this strong increase in the magnetic pressure in 
front of the cloud turns out to be a result rather independent of the initial cloud geometry. In 
fact, Maxwell stresses will tend to squeezed the cloud into a cylinder- like structure with the axis 
aligned perpedicular to the plane defined by the initial field and the direction of motion. In this 
respect, low resolution simulations started with a cylinder as the initial cloud shape, revealed a 
qualitatively similar behavior in the magnetic energy evolution, although an even stronger increase 
in the magnetic pressure occured at the cloud nose. Overall these results compare well with the 
previous 2-D simulations of Miniati et al. (1999b). While the increase of magnetic energy in the 3-D 
simulation is slightly lower than for the 2-D case, the values are comparable in order of magnitude. 



3.2.3. Timescales 

In the evolution of the cloud several dynamical timescales can be identified. As the magnetic 
pressure in front of the cloud increases, the flow progressively tends to be dominated by Maxwell 
stresses and its global behavior radically changes from the pure hydro dynamical case. Typically, 
this occurs when the magnetic pressure becomes comparable to the ram pressure of the gas. Gregori 
et al. (1999) have shown that the onset of this transition is characterized by magnetically enhanced 
R-T modes that ultimatively disrupt the entire cloud. In Fig. 6 the ratio of the maximum magnetic 
pressure, ps, along the symmetry axis (y = z = 0) with respect to the ram pressure of the intercloud 
fluid, PiUc, is plotted as a function of time for the two values of (3. At the beginning of the simulation 
PBo/(Pi u c) ~ 0.067 for (3 = 4 and p bo /(piU^) ~ 0.003 for (3 = 100. Thus, the ram pressure is initially 
dominant. However, as shown in Fig. 6 for the (3 = 4 simulation, at about one crushing time the 
magnetic pressure maximum becomes comparable and even larger than the initial ram pressure. Its 
subsequent decrease is connected to the rapid development of the R-T instability; in fact the field 
lines move away from the symmetry axis, towards the newly developing indentations (see Fig. 2). 
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Our analysis of the field inside such structures reveals that the magnetic pressure keeps growing 
there and at t ~ 2r cr reaches a ratio of 1.66 with respect to the ram pressure. For the (3 = 100 case 
we note that, despite a remarkable increase with respect to its initial value, the magnetic pressure 
maximum remains well below the ram pressure limit throughout the simulation time. Based on 
the model developed in the previous sections, after the initial transient, from equation (3-9) we 
estimate 

^ikir) ' (3 ~ 12) 
where Ma is the initial Alfenic Mach number. Using m ~ 1, we estimate that the magnetic pressure 
becomes comparable to the initial ram pressure for t > (M]/2)»r s ~ (M\/2)t s , or t > 0.28r cr 
((3 = 4) and t > 7.03r cr ((3 = 100). This is indeed confirmed from Fig. 6. For the (3 = 4 simulation 
at t > 0.3r cr the magnetic pressure maximum is already comparable with the ram pressure. On the 
other hand, in the weak field case, the magnetic pressure maximun remains much below the initial 
ram pressure even at simulation end. These results roughly confirm our findings for the behavior 
of the magnetic field at the cloud nose. 

Another significant timescale can be defined in terms of the energy exchange between kinetic 
and magnetic energy, since the stretching of the field lines acts as a conversion mechanism in which 
the cloud kinetic energy is transformed into magnetic energy. By analogy with Miniati et al. (1999b) 
we can estimate the time required for the magnetic energy to equal the initial kinetic energy of the 
cloud by solving the equation 

1 2 „ „ „ „ ( t~ 



2 p c u 2 c V c = E B - E B0 ~ PB0 V C \-j . (3-13) 

We easily get the characteristic timescale 

f X M A \ 1/m f X M A \ 

where the last term has been obtained by setting m = 1, as previously observed. In this respect, T ma 
must be intended as a characteristic braking time, since it sets an upper limit for the cloud motion 
across the ISM. We should note that in absence of magnetic field, the characteristic time required to 
stop the cloud is given by Td = X T s (e.<?., Jones et al. 1994). It is then evident that the cloud motion 
is magnetically dominated only if r ma < r^. By comparison with (3-14), we see that this is the 
case if Ma < 4. This indicates that for (3 = 4 (Ma = 2.74) the magnetic field plays a dynamically 
important role from the beginning. On the other hand, in the (3 = 100 case, Ma = 13.69 and the 
hydrodynamic drag is the relevant timescale in determining the cloud evolution; thus in this case 
magnetic effects are likely not to be important. Mouschovias & Paleologou (1979) and Elmegreen 
(1981) have calculated the characteristic time required by the magnetic field to stop the cloud using 
a different approach based on the radiated energy through Alfven waves. Typically, for transverse 
field geometry and an infinitely long cylindrical cloud, they estimate that the timescale for magnetic 
braking is of the order of \ \Mat s which, apart from a numerical factor, 2, is equivalent to our 
result. 



- 11 - 



The timescale (3-14) can be recast in a dimensional form: 



T~ma — 



7.04 x lO 10 ^ 1 / 2 




(3-15) 



The characteristic time for radiative cooling is given by (Spitzer 1978) 



3k B T 



~ 1.6 x 10 i4 s 



(3-16) 



2 rijA 



where rii is the number density in the intercloud medium and A is the interstellar cooling function. 
For temperatures of the interstellar gas in the order of 10 4 K and densities rii ~ 0.1 cm" 3 , we have 
njA ~ 1.3 x 10~ 26 erg/s (Ferrara & Field 1994) to give the numerical value in equation (3-16). We 
clearly see that for the values of the parameters x an d (3 used in the simulations, the evolution of 
the system remains always strictly adiabatic (that is, T ma < t coo \) only for for clouds with radius 
R c < 0.1 pc. In the case of larger clouds, radiative effects have still a small influence on the cloud 
dynamics if the cooling time is larger than the characteristic time for re-expansion of the postshock 
interstellar material as it flows around the cloud. Such time is in the order of R c /c s (Miniati et 
al. 1997). From the condition T ma > R c /c s we then obtain that cooling remains negligible for clouds 
radii R c < 10 pc. 



When the magnetized fluid flows into the cloud wake, the field lines frozen in it produce 
elongated structures of strong field concentration surrounded by a thin vortex and current sheet. 
These filamentary structures are clearly visible in Figs. 7 and 8. Following Mac Low et al. (1994) 
we identify these regions as "flux-ropes" , where individual flux tubes are organized in a coherent 
pattern. In the simulations presented here, however, we see no signs of twist in them. The complete 
3-D structure is illustrated in Fig. 9 for t = 0.94r cr . In both the the strong and the weak 
field simulations, the flux ropes converge toward the symmetry axis where the thermal pressure 
is low (Mac Low et al. 1994). There are then strong similarities between our results and the 2- 
D simulations of Mac Low et al. (1994) and Jones et al. (1996), which show the same flux rope 
structures for the transverse field geometry. From Fig. 7 we can see that in the weak field {[3 = 100) 
simulation the flux rope remains very close to the symmetry axis for almost its entire length. Inside 
the flux rope the field strength at this time (t = 0.94r cr ) is comparable to that at the leading edge of 
the cloud. In the cloud wake the dominant enhancement of the field strength is probably produced 
by the compression of the rope on the symmetry axis. Such compression, in turn, is due to the 
large gradient of thermal pressure produced there by the cloud motion. We should also note that 
some additional field amplification may occur as the fluid crosses the tail shocks. Conversely, in 
the strong field, (3 = 4, simulation (Fig. 8), the flux rope is kept close to the symmetry axis only in 
a small region just at the cloud rear (a few cloud diameters in size). Farther away it opens up in 
"wings" that tend asymptotically to align with the initial unperturbed field direction. Finally, the 
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magnetic field value (at t = 0.94r cr , = 4) in these ropes is on average a factor ~ 2 larger than 
the background field. 

As we can see in Fig. 10, an important feature of the flux ropes is that the plasma in them is 
dynamically tied to the cloud; that is, the plasma there moves with the cloud. This is a very re- 
markable characteristic because the flux ropes, formed by magnetic field lines stretched around the 
shape of the cloud, do not carry cloud material per se (demonstrated through the tracer variable C). 
The separation of the flux ropes from the external flow occurs through a tangential discontinuity, 
with the magnetic field parallel to the interface and no mass flow across it. Such a discontinuity is 
characterized by a constant total pressure (magnetic plus thermal) across the separation interface. 
However, the higher magnetic pressure inside the rope is balanced by a considerably lower density 
there, compared to the outside fluid material. As already noted, there is little or no twist generated 
by magnetic or current helicity in the flux ropes. They are sheet-like. Similar sheet-like morpholo- 
gies were noted in strong magnetic structures formed in high resolution 3-D MHD Kelvin-Helmholtz 
instability simulations (Ryu et al. 2000). In that case it was clear that the sheet-like morphology 
was a consequence of the strong magnetic field and not a resolution effect. This may be true here 
as well, but our resolution is not sufficient to assess that. 



3.4. Magnetic Reconnection 

The ideal MHD approximation breaks down in regions where resistive effects or diffusion dom- 
inate the field evolution. Even if the scale length of these regions is very small, their development 
may change the global topology of the field lines (Song & Lysak 2000). The most important of 
such events is magnetic reconnection. This has been observed in the heliosphere and in various 
2-D and 3-D numerical simulations (e.g., Biskamp 1994, Miniati et al. 1999b, Antiochos & De Vore 
1999, Ryu et al. 2000). Miniati et al. (1999b), using the 2-D version of the MHD code applied here, 
showed clear 2-D reconnection in the wake of a cloud, apparently due to tearing mode instabilities. 
Our new numerical simulations confirmed the previous results for clouds and showed that, indeed, 
reconnection takes place in the wake of the cloud. In fact from Figs. 7 and 8, it is clear that on 
the symmetry plane, y = 0, a pair of thin current sheets form and that field annihilation occurs 
there. Here, the reconnection process seems steady, however, that can be understood by examining 
the properties of the reconnection region. The two figures show that the ratio L/5 between the 
length and the full thickness of the current sheet is ~ 15-25 for both the strong and weak field 
simulations. The current sheet structure, for the (3 = 4 case, reveals strong similarities with the 
classic 2-D Sweet-Parker current sheet model (e.g., Biskamp 1993). In this respect, the recent 2-D 
resistive MHD simulations of Uzdensky & Kulsrud (1999) also seem to confirm that reconnection 
tends spontaneously to evolve toward a Sweet-Parker steady state. Additional analogies with other 
reconnection models, e.g., Syrovatskii's solutions (Biskamp 1993), can also be found in the branch- 
ing of the "separatrix" on the right side of the sheet. In the presence of a reconnection layer the 
issue regarding its stability immediately is raised. Stability against tearing modes is expected if 
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S ~ {L/5) 2 < 10 4 (Biskamp 1993), where 5 is the Lundquist number of the current sheet. This is 
indeed in agreement with our present simulations, which show no development of such an instability. 
That contrasts as expected with our previous 2-D simulations where the tearing mode instability 
was, in fact, observed in the analogous regions Miniati et al. (1999b). In order to understand the 
difference between the two cases, we recall that the Lundquist number is actually defined as 

S = ^, (3-17) 

where va is the Alfven velocity of the magnetic field on the top and bottom of the current sheet and 
T) is the fluid resistivity. In the Sweet-Parker reconnection model S takes the form S = (L/5) 2 . In 
simulations of ideal MHD 5 is limited by numerical discretization and its value is ~ 2-3 Ax, where 
Ax is the size of a numerical zone. The length of the reconnection region in the cloud wake will 
generally be several cloud radii, typically 5-10. This implies that the values for L/5 are numerically 
limited: S < (2-3 x N c ) 2 , where N c is the number of numerical zones per cloud radius. In our 
current 3-D simulations, N c = 26, so S < 10 4 , and the layer should be stable to tearing modes as, 
indeed, it is. On the other hand, our 2-D simulations were done at roughly twice the resolution 
across the clouds, with N c = 50, so the current sheet in the tail was characterized by S > 10 4 . It 
should have been, and was, tearing mode unstable. 



4. Summary & Concluding Remarks 

In this paper we have presented the results of a series of 3-D MHD numerical simulations of 
cloud-wind interactions. We have considered an initial spherical cloud that moves transverse to 
the magnetic field, with two different cases for its initial strength; namely (3 = p/pb = 4 and 
(3 = 100. Both the weak (f3 = 100) and strong field ((3 = 4) simulations showed a qualitatively 
comparable behavior with a substantial enhancement of the magnetic pressure at the leading edge 
of the cloud. This confirms and extends previous 2-D results (Jones et al. 1996, and Miniati et 
al. 1999b). A new, detailed analysis of the field stretching that occurs in front of the cloud is 
presented here. Qualitative agreement with the numerical simulations is obtained. In particular, 
it is shown that magnetic effects tend to be important only if the initial Alfvenic Mach number is 
sufficiently small. In this case the cloud dynamical evolution (acceleration and disruption) occurs 
on shorter timescales and a substantial rapid conversion of kinetic into magnetic energy takes place. 
Thus, at t < r cr , the magnetic pressure becomes comparable with or greater than the ram pressure. 
On the other hand, if the initial Alfvenic Mach number is large, then the evolution of the cloud is 
still dominated by the hydrodynamic drag and conversion of kinetic into magnetic energy occurs 
at a much smaller rate. 

For typical values of the magnetic field in the ISM, the propagation of a cloud through a 
magnetized medium should produce inhomogeneities on a scale comparable to the cloud size. As 
already pointed out in Miniati et al. (1999b), this could help explain observed fluctuations in the 
galactic magnetic field (Heiles 1989, Meyers et al. 1995). Observations of Gloeckler et al. (1997) 
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in the local interstellar cloud also seem to suggest the possibility of an inhomogeneous distribution 
of the field intensity. Moreover, we illustrate a complex series of topological modifications that 
take place in the cloud wake and characterize the dynamics of the flow there. In the tail of 
the cloud the field lines aggregate in a coherent pattern to form long, sheet-like flux ropes. The 
plasma entrained there moves approximately with the cloud. Additionally, along the symmetry 
axis, magnetic reconnection takes place. Evolutionary details such as the formation of vortical 
structures at the cloud interface, as well as the development of instabilities of the reconnection 
current sheet were prevented by the limited numerical resolution. Despite these limitations, the 
simulations reported here give a significant picture of the fundamental processes governing cloud- 
wind interaction in a magnetized ISM, setting important constraints for further observational and 
theoretical studies. In this respect, our results may give support to recent models of cloud-wind 
interactions such as those proposed to explain the formation of nonthermal filaments in the galactic 
center (Shore & LaRosa 1999). Investigations of planetary nebulae interacting with a magnetized 
fluid may also take advantage of the results presented here (see e.g., Soker & Dgani 1997, Dgani & 
Soker 1998). 
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FIGURE CAPTIONS 

Fig. 1. — Schematic of the numerical setup. 

Fig. 2. — Slice of log(p) on the plane y = at t = 0.94r cr (left column) and t = 2.25r cr (right 
column). Top: (3 = oo, center: (3 = 100, bottom: (3 = 4. 

Fig. 3. — Relative cloud momenta in the y and z directions (see text). 

Fig. 4. — Top: magnetic energy increment in the computational volume normalized with respect 
to the initial magnetic energy inside the cloud (AEb/E B0 ). Bottom: magnetic energy increment 
in the computational volume normalized with respect to the initial kinetic energy of the cloud 
(AEb/Ekq). The high resolution (hr) simulation corresponds to case 3 in Table 1 and the low 
resolution (lr) simulation corresponds to case 4 in Table 1. 

Fig. 5. — Total kinetic energy in the computational volume normalized with respect to its initial 
value (Ek/Ekq). The high resolution (hr) simulation corresponds to case 3 in Table 1 and the low 
resolution (lr) simulation corresponds to case 4 in Table 1. 

Fig. 6. — Maximum value of the magnetic pressure on the symmetry axis (y = z = 0) with respect 
to its initial ram pressure of the fluid. 

Fig. 7. — Slice in the plane z = for the (3 = 100 simulation at t = 0.94r cr . log(p) top left, 
log(i? 2 /2) top right, log(w 2 + 1) bottom left, log(j 2 + 1) bottom right. Here, u = V x u is the 
vorticity and j = V x B is the current density. The B x and B y components of the magnetic field 
are represented by the arrows. 

Fig. 8. — Slice in the plane z = for the (3 = 4 simulation at t = 0.94r cr . log(p) top left, log(-B 2 /2) 
top right, log(c<j 2 + 1) bottom left, log(j 2 + 1) bottom right. Here, u) = V x u is the vorticity and 
j = V x B is the current density. The B x and B y components of the magnetic field are represented 
by the arrows. 

Fig. 9. — Left column: volume rendering of the magnetic pressure (log scale) for the (3 = 100 
simulation (top) and (3 = 4 (bottom). Right column: volume rendering of the cloud density (log 
scale) for the (3 = 100 simulation (top) and (3 = 4 (bottom). Time is expressed in units of r cr . 

Fig. 10. — Slice of log(p) on the plane z = at t = 0.94r cr (left column) and t = 2.25r cr (right 
column). Top: (3 = oo, center: (3 = 100, bottom: (3 = 4. The u x and u y components of the velocity 
field are represented by the arrows. 
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Table 1. Summary of 3-D MHD cloud simulations. 



Case 


Symmetry 
plane(s) a 


Resolution 13 


Initial field 
direction 


Shape 

(Ixi lyi lz) 


M 


X 


p 


M A 


1 


x-y, y-z 


26 


(0,1,0) 


(1,1,1) 


1.5 


100 


oo 


CO 


2 


x-y, y-z 


26 


(0,1,0) 


(1,1,1) 


1.5 


100 


100 


13.7 


3 


x-y, y-z 


26 


(0,1,0) 


(1,1,1) 


1.5 


100 


4 


2.74 


4 


x-y, y-z 


16 


(0,1,0) 


(1,1,1) 


1.5 


100 


4 


2.74 


5 


none 


16 


(A. JL J-) 


(1,1,2) 


3.0 


100 


10 


8.66 


6 
7 
8 
9 


none 
none 
none 
none 


16 
16 
16 
16 


(A. JL J_) 


(1,1,2) 
(1,1,2) 
(1,1,2) 

(1,1,2) 


3.0 
3.0 
3.0 
3.0 


10 
30 
100 
100 


10 
10 

2 
4 


8.66 
8.66 
3.87 
5.48 



a plane(s) across which symmetric boundary conditions are employed in the numerical 
computations. 

b in number of zones per cloud radius, R c . 

c initial cloud axis in units of R c . (1,1,1) corresponds to a sphere and (1,1,2) corre- 
sponds to a cylinder with the axis in the z direction. 
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